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Abstract 

One-hundred-and-three isolates of Bacteroides ovatus, B. thetaiotaomicron, and 
5. xylanisolvens were recovered from cow, goat, human, and pig fecal enrich- 
ments with cellulose or xylan/pectin. Isolates were compared using 16S rRNA 
gene sequencing, repetitive sequence-based polymerase chain reaction (rep- 
PCR), and phenotypic microarrays. Analysis of 16S rRNA gene sequences 
revealed high sequence identity in these Bacteroides; with distinct phylogenetic 
groupings by bacterial species but not host origin. Phenotypic microarray analy- 
sis demonstrated these Bacteroides shared the ability to utilize many of the same 
carbon substrates, without differences due to species or host origin, indicative 
of their broad carbohydrate fermentation abilities. Limited nitrogen substrates 
were utilized; in addition to ammonia, guanine, and xanthine, purine deriva- 
tives were utilized by most isolates followed by a few amino sugars. Only rep- 
PCR analysis demonstrated host-specific patterns, indicating that genomic 
changes due to coevolution with host did not occur by mutation in the 16S 
rRNA gene or by a gain or loss of carbohydrate utilization genes within these 
Bacteroides. This is the first report to indicate that host-associated genomic dif- 
ferences are outside of 16S rRNA gene and carbohydrate utilization genes and 
suggest conservation of specific bacterial species with the same functionality 
across mammalian hosts for this Bacteroidetes clade. 



* Mention of a trade name, proprietary 
product, or specific equipment does not 
constitute a guarantee or warranty by the 
USDA and does not imply approval to the 
exclusion of other products that may be 
suitable. 



Introduction 

The gastrointestinal tracts (GIT) of vertebrates are among 
some of the most densely populated environments for 
microbes, with cell counts of -lO^^-lO^^ cells per mL 
(Ley et al 2006; Thomas et al 2011). These microbial 
populations play an important role in the development of 
the immune system, health, and nutrition of the verte- 
brate host (Dethlefsen et al. 2007). The advent of 16S 
rRNA gene sequence phylogeny has allowed for analysis 
of whole microbial communities. Intestinal bacterial pop- 



ulations show a low diversity at the phyla level, domi- 
nated by Firmicutes and Bacteroidetes; other phyla 
represented include Protobacteria, Actinobacteria, Fusobac- 
teria, Verrucomicrobia, Fibrobacteres, and Spyrochaeates 
(Eckburg et al. 2006; Ley et al. 2008a). At the species and 
strain levels, the amount of diversity increases dramati- 
cally which has been partially attributed to host selective 
pressure and coevolution (Backhed et al. 2005; Ley et al. 
2006). 

The intestinal bacteria and their hosts live together in 
mutualistic relationships. The microbes that reside in the 
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GIT have the capacity to produce enzymes that degrade 
dietary polysaccharides from plant cell walls that the 
host's enzymes cannot (Salyers 1984; Flint et al. 2008). 
The fermentation of these polysaccharides produces end 
products, such as short-chain fatty acids (SCFA), that the 
host uses as nutritional sources (Martens et al. 2011). The 
amount of energy supplied by the SCFA produced can be 
as much as 70% in ruminants, 20-30% in many omnivo- 
rous animals, and 10% in humans (Bergman 1990). 

The intricate mutualistic relationship of the microbiota 
in the distal gut and the host results in an environment 
that allows for the coevolution of both microbiota and 
host. This coevolution, defined as the reciprocal adapta- 
tions of each lineage in response to the other, has resulted 
in differences in microbial composition among host spe- 
cies and is helpful in understanding the relationships 
between host and microbiota (Zaneveld et al. 2009). 
Comparison of microbial community 16S rRNA gene 
sequences from fecal samples across a wide range of host 
species revealed that microbial communities codiversified 
with their host (Ley et al. 2008b); suggesting evolutionary 
interaction between a host and their microbiota (Ley 
et al. 2008b). Coevolution could also be investigated 
using a single bacterial species or group that is found in 
multiple host species. These bacteria could be compared 
directly for shared ancestry; one would expect to find 
subpopulations specific to each host (Oh et al 2010). 

Eight week continuous culture enrichments of feces 
with cellulose and xylan/pectin (2/1) were used to isolate 
bacteria from the polysaccharide degrading microbial 
populations using feces from four hosts: cows {n = 4), 
goats (n = 4), humans (« = 4), and pigs (n = 6) (Ziemer 
2013). A total of 1650 bacteria were isolated and identi- 
fied using 16S rRNA gene sequencing; phyla distribu- 
tions were similar across hosts; averaging 45.1% 
Firmicutes, 32.7% Bacteroidetes, 12.4% Proteobacteria, 
2.6% Actinobacteria, 5.5% Fusobacteria, and 0.7% Syn- 
ergistestes. These results are similar to those reported by 
Ley et al. (2008a) for feces, in 60 different mammalian 
species, Firmicutes and Bacteroidetes made up 82% of 
the nearly 20,000 classified sequences with 14.2% in 
the Proteobacteria, Fusobacteria, and Actinobacteria 
phyla. Nearly 75% (402 of 540) of Bacteroidetes isolates 
were identified as B. ovatus, B. thetaiotaotnicron, or B. 
xylanisolvens, making this one of the predominant 
groups. Multiple isolates with slightly different 16S rRNA 
gene sequences were isolated from enrichments from 
each fecal donor. 

Subgroups of Bacteroides species have been found 
within different hosts suggesting that they coevolved 
together (Backhed et al. 2005). The gastrointestinal Bacte- 
roides are most often transmitted from parent to offspring 
(Ley et al. 2006). Bacteroides species are anaerobic, non- 
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motile, gram-negative rods and are one of the most 
numerous bacteria found in the colon of many different 
animal species (Thomas et al. 2011). Bacteroides have an 
unprecedented ability to degrade both plant poly- 
saccharides and host mucins. A predominant group of 
Bacteroides species, including B. ovatus, B. thetaiotaotnicron, 
and B. xylanisolvens, from multiple hosts (pig, cow, goat, 
and human) were characterized using 16S rRNA gene 
phylogenetic analysis, repetitive sequence-based PCR 
(rep-PCR) genomic banding patterns, and phenotypic 
microarrays of carbon and nitrogen metabolism to deter- 
mine if differences due to host origin could be detected. 

Experimental Procedures 
Isolations 

Bacteroides isolates were obtained from enrichments of 
fecal samples taken from four cows, four humans, four 
goats, and six pigs. Freshly voided feces (at least 500 g) 
were collected from individuals and returned to the lab 
for processing within 30 min. Continuous culture enrich- 
ments with cellulose and xylan-pectin mixture were per- 
formed for 8 weeks with each fecal sample as described in 
Ziemer (2013). Host type is designated in the isolate 
names as follows: NLAE-zl-Cx for isolates from cow, 
NLAE-zl-Gx for isolates from goat, NLAE-zl-Px for iso- 
lates from pig or NLAE-zl-Hx for isolates from human 
fecal enrichment cultures (see Table SI for NCBI Acces- 
sion Numbers). The isolates were fairly evenly distributed 
across individuals averaging five per cow, eight per goat, 
seven per human, and six per pig. Type strains, B. ovatus 
ATCC 8483 (X83952) (ATCC, Manassas, VA), B. thetai- 
otaomicron VPI 5482 (L16489) (ATCC), and B. xylanisolvens 
XBIA (AM230650) (DSMZ, Braunschweig, Germany) 
were included in analyses. A total of 85 cow, 126 goat, 77 
human, and 114 pig Bacteroides, within B. ovatus, B. 
thetaiotaomicron, and B. xylanisolvens species clade were 
isolated in pure culture and DNA was extracted using 
DNeasy kit (QIAGEN Sciences, Germantown, MD). These 
isolates were analyzed by 16S rRNA gene sequencing and 
repetitive sequence-based PCR (rep-PCR) as described 
below. AU of the 16S rRNA gene sequences were within 
95% similarity regardless of Bacteroides species with no 
host-specific clustering. Isolates with rep-PCR banding 
patterns that were unique within a Bacteroides species 
were selected for further analysis. When rep-PCR patterns 
were identical among isolates within a Bacteroides species 
from a single fecal donor, the isolates were considered 
clonal and one isolate was selected as representative. 
When isolates from the same host type but different fecal 
donors demonstrated similar rep-PCR patterns then one 
isolate representing each fecal donor was chosen. The 
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final number of isolates for analysis was 103 total: 17 B. 
ovatus, 32 B. thetaiotaomicron, and 54 B. xylanisolvens (21 
cow, 31 goat, 29 human, and 22 pig isolates). 

16S rRNA gene sequencing 

16S rRNA gene sequencing of isolates used the 27f and 
1492r primers (Lane 1991) in a 50 /,(L reaction. The reac- 
tions contained 5 jiL 10 x PGR buffer with 15 mmol/L 
MgClj, 0.5 ^(L of 10 mmol/L dNTPs, 1 ^L each of 
25 pmol/L primer, 0.125 /iL of 5 U//iL Taq polymerase 
(Invitrogen, Carlsbad, CA), 40 /iL of nuclease-free water, 
and -100 ng of sample DNA. A MJ Research PTC-220 
thermocycler (MJ Research, St. Bruno, Quebec, Canada) 
was used with the following program: 95°C for 5 min 
and then for 30 cycles at 95°C for 45 sec, 47°C for 1 min, 
72°C for 1 min, followed by a 72°C for 7 min final exten- 
sion. Sequencing was done on a 3730x1 DNA Analyzer 
(Applied Biosystems, Carlsbad, CA) by the Iowa State 
University DNA Sequencing and Synthesis Facility (Ames, 
lA). 

Consensus sequences were generated using VNTI 11.1 
software (Invitrogen), and those >1200 bp in length were 
analyzed using Ribosomal Database Project II (Cole et al. 
2009) Sequence Match for identification. AH sequences 
have been deposited in GenBank (Accession numbers 
presented Table SI). For this experiment, bacteria that 
identified at >90% B. ovatus, B. thetaiotaomicron, and 
B. xylanisolvens were selected for further investigation. 
Sequences identified with Bacteroides strain designations 
were also included if homologous to the targeted clade. 
Sequence similarity was analyzed using Bionumerics 
(v. 6.5 software. Applied Maths, Austin, TX) and a 
dendogram created using standard pairwise alignment 
and unweighted pair group method with arithmetic mean 
(UPGMA) clustering. 

Repetitive sequence-based PCR 

The bacterial isolates within the B. ovatus - B. thetaiota- 
omicron — B. xylanisolvens clade were analyzed using rep- 
PCR with BOXAIR primer and ERICIR and ERIC2 prim- 
ers (Versalovic et al. 1991). The rep-PCR reactions for 
the BOX primer contained 3.35 /(L 10 x PCR buffer with 
15 mmol/L MgCh, 0.5 /iL of 10 mmol/L dNTPs, 1 /iL of 
25 pmol/L BOXAIR primer, 0.125 /xL 5 U//(L Taq poly- 
merase (Invitrogen), and -100 ng of sample DNA in 
25 iiL volume. The reactions for ERICIR and ERIC2 were 
the same except that 5 /iL 10 x PCR buffer with 
15 mmol/L MgClz and 1 /iL each of 25 pmol/L ERICIR 
and ERIC 2 primers, were used per 25 iiL reaction. All 
reactions were run on an MJ Research PTC-220 thermo- 
cycler (MJ Research). The BOX and ERIC reaction condi- 



tions were 95°C for 5 min and then for 30 cycles at 95°C 
for 45 sec, 50°C and 46°C, respectively, for 1 min, 72°C 
for 1 min, followed by 72°C for 7 min final extension. 
The PCR amplicons were separated by electrophoresis on 
a 1.5% Agarose (Amresco, Solon, OH) gel in Ix TBE 
buffer for 6 h at 100 V with 100 bp molecular weight 
marker (Amresco). Gels were stained with ethidium bro- 
mide and pictured under ultraviolet light using the Kodak 
Image Station 4000MM Pro (Molecular Imaging Systems, 
Garestream Health Inc., Rochester, NY). Analysis of the 
gel banding patterns was done using BioNumerics 
(Applied Maths) and dendograms were generated using 
Pearson correlation similarity coefficient with optimiza- 
tion of 1% and the UPGMA clustering method. 

Phenotypic comparison 

The 103 isolates were then screened for their ability to 
metabolize 190 carbon and 95 nitrogen sources using 
phenotypic microarrays (PMl and PM2A for carbon and 
PM3B for nitrogen) (Biolog, Hayward, CA), according to 
manufacturer's instructions without modifications (Mapp- 
ley et al. 2012). Microarray plates were inoculated and 
incubated anaerobically for 24 h then read using a Bio- 
log MicroStation reader. Growth was measured using an 
optical density at 590 nm for PMl and PM2A an optical 
density at 750 nm for PM3B. The following carbon 
sources were not analyzed due to abiotic reactions 
under anaerobic experimental conditions: arabinose, 
dihyderoxyacetone, D-gluosamine, 2-hydroxybonzoic 
acid, 5-keto-D-gluconate, L-lyxose, palatinose, D-ribose, 
3-deoxy-D-ribose, sorbic acid, D-tagatose, and D-xylose 
(Borglin et al. 2009; Mappley et al. 2012). The values 
were zeroed to account for background noise and any 
results <0.1 O.D. (>2 times background O.D.) were 
considered to have no growth. In order to assess repeat- 
ability of the phenotypic arrays a subset of isolates were 
analyzed on two to five plates. The two carbon utiliza- 
tion arrays were combined for analysis of utilization 
patterns. Analysis of the data was done using BioNumer- 
ics (Applied Maths) for carbon and nitrogen substrate 
sets. The dendograms were generated using the Pearson 
correlation similarity coefficient and UPGMA clustering 
method. 

Results 

16S rRNA phylogeny 

Analysis of 16S rRNA gene sequences for B. thetaiotaomi- 
cron, B. ovatus, and B. xylanisolvens isolates and type 
strains (ATCC, Manassas, VA; DSMZ, Braunschweig, 
Germany) demonstrate all were within >97% similarity to 
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B. ovatus: NLAE-ZI-H251, NLAE-zl-H255, NLAE-zl-H356, S. ovatus ATCC 8483, NLAE- 
ZI-H304, NLAE-ZI-H361, NLAE-zl-H163, NLAE-zl-H426, NLAE-zl-H307, NLAE-zl-H545, 
NLAE-ZI-H313, NLAE-zl-H541, NLAE-zl-H422 



I<;3 B. ovatus: NLAE-zl-Cll, NLAE-zl-C34 



B. xylanisolvens: NLAE-zl-C29, NLAE-zl-H194, NLAE-zl-H465, NLAE-zl-G481, NLAE-zl- 
G445, NLAE-ZI-G435, NLAE-zl-G405, NLAE-zl-G515, NLAE-zl-P218, NLAE-zl-P718, 
NLAE-ZI-P371, NLAE-ZI-P402, NLAE-zl-P225, NLAE-zl-H40, NLAE-zl-G157, NLAE-zl- 
G310, NLAE-ZI-P299, NLAE-zl-G421, NLAE-zl-P352, NLAE-zl-P393, NLAE-zl-P727, 
NLAE-ZI-P735, NLAE-zl-P732, NLAE-zl-C139, NLAE-zl-C195, NLAE-zl-C233, NLAE-zl- 
C315, NLAE-ZI-C339, NLAE-zl-G346, NLAE-zl-C182, NLAE-zl-G397, NLAE-zl-GllO, 
NLAE-ZI-G353, NLAE-zl-C178, NLAE-zl-P295, NLAE-zl-C257, NLAE-zl-G37, NLAE-zl- 
G39, NLAE-ZI-G200, NLAE-zl-G237, NU\E-zl-G275, NLAE-zl-G193, NLAE-zl-G194, 
NLAE-ZI-P349, NLAE-zl-G109, NLAE-zl-C367, NLAE-zl-G408, NLAE-zl-G2, 
NLAE-ZI-G44, NLAE-zl-P887, NLAE-zl-G228, B. xylanisolvens XBIA, NLAE-zl- 
P412, NLAE-ZI-G289 



NLAE-ZI-G312 xylanisolvens 
e. thetaiotaomicron: NLAE-zl-P696, NLAE-zl-P737, NLAE-zl-P745 



S. thetaiotaomicron: NLAE-zl-P174, NLAE-zl-P750, NLAE-zl-C523, NLAE-zl-C13, 
NLAE-ZI-C516, NLAE-zl-C504, NLAE-zl-C557, NLAE-zl-G295, NLAE-zl-G303, NLAE-zl- 
G288, NLAE-ZI-H492, NLAE-zl-H463, NLAE-zl-P32, NLAE-zl-H353, NLAE-zl-H39, 
NLAE-ZI-H373, NLAE-zl-H421, NLAE-zl-P699, NLAE-zl-G234, NLAE-zl-H4S6, NLAE-zl- 
H23, NLAE-ZI-H536, NLAE-zl-H409, NLAE-zl-H334, NLAE-zl-G493, NLAE-zl-C15, 
NLAE-ZI-H207, S. tlietaiotaomicron\/P\ 5482, NLAE-zl-C484 



NLAE-ZI-C425 
NLAE-ZI-C501 
NLAE-ZI-H59 
NLAE-ZI-H73 



thetaiotaomicron 
ovatus 
ovatus 
ovatus 



Figure 1. Cluster analysis of 16S rRNA genes using UPGMA, dendogram with similarity values (%). Bar at top is similarity continuum. 
Differentiation among host by isolate name: cow = NLAE-zl-Cx, goat = NLAE-zl-Gx, human = NLAE-zl-Hx, and pig = NLAE-zl-Px. 
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each other (Fig. 1). There were three exceptions, all B. 
ovatus strains, one cow isolate (NLAE-zl-C501) at 96.1%, 
and two human isolates (NLAE-zl-H59 and NLAE- 
zI-H73) at 94.9% similarity with the other isolates. The 
isolates formed distinct phylogenetic groupings by bacte- 
rial species with 17 isolates grouped with B. ovatus, 32 
with B. thetaiotaomicron, and the remaining 54 with 
B. xylanisolvens type strains. Bacteroides ovatus strains 
were isolated from human (14) and cow (3) fecal enrich- 
ments. These isolates distinctly grouped by host origin 
with 98% sequence similarity. The B. thetaiotaomicron 
isolates (seven pig, eight cow, five goat, and 12 human) 
had >98.6% sequence similarity with no differentiation 
due to host origin. Similarly, the B. xylanisolvens isolates 
(15 pig, 10 cow, 26 goat, and three human) had >98.9% 
sequence similarity, again, no effects of host origin. 

Analysis of the rep-PCR banding demonstrated dis- 
tinct patterns with both the BOX and ERIC primers by 
Bacteroides species and host origin, results are presented 
by Bacteroides species (Figs. 2, 3). Within B. ovatus, iso- 
late banding patterns from BOX rep-PCR formed three 
groups of human isolates with >90% similarity except 
for isolate NLAE-zl-H541 which had 89.4% similarity 
with the type strain (S. ovatus ATCC 8483) (Fig. 2A). 
The BOX rep-PCR banding patterns of the two human 
(NLAE-Z1-H59 and NLAE-zl-H73) and one cow (NLAE- 
zl-C501) isolates with less than 97% 16S rRNA gene 
sequence similarity were distinct from the other B. ovatus 
isolates. The ERIC rep-PCR produced four groups 
of B. ovatus human isolates with >90% similarity 
(Fig. 3A). The banding patterns of NLAE-zl-H59, 
NLAE-zl-H73, and NLAE-zl-C501 were more similar to 
the rest of the isolates using the ERIC primers in con- 
trast to results using BOX rep-PCR. In fact the cow iso- 
late, NLAE-zl-C501, has the highest banding pattern 
similarity to the type strain B. ovatus ATCC 8483 at 
82.8%. 

Analysis of banding patterns generated using the BOX 
rep-PCR for B. thetaiotaomicron isolates at >90% similar- 
ity formed one group of goat isolates, four groups of 
human isolates, one group of pig isolates, and one group 
of a mix of pig and cow isolates. Within this mixed 
group, the type strain B. thetaiotaomicron VPI 5482 is also 
found (Fig. 2B). The ERIC primers produced one cow, 
three human, and three mixed groups with >90%. Two of 
the mixed groups contain human and cow isolates, where 
the third group contains goat, human, and pig isolates. 
The B. thetaiotaomicron VPI 5482 type strain and NLAE- 
zl-H207 did not group with any of the other isolates 
(Fig. 3B). 

Repetitive sequence-based PCR of B. xylanisolvens using 
the BOX primer resulted in very few bands, ranging in 
number between 1 and 4, thus was not used to assess B. 



xylanisolvens genetic diversity. Using a >90% similarity 
cutoff, ERIC rep-PCR banding pattern analysis resulted in 
two goat, two pig, and one mixed, cow and goat, isolate 
groupings. The type strain B. xylanisolvens XBIA was 
most similar to a human isolate, NLAE-zl-H465, at 
91.9%. Four isolates (NLAE-zl-H194, NLAE-zl-H40, 
NLAE-Z1-P412, and NLAE-zl-G481) had distinct banding 
patterns compared to the other B. xylanisolvens isolates 
(Fig. 3C). 

Substrate utilization 

Analysis of isolate carbon and nitrogen utilization pat- 
terns produced with the phenotypic microarrays did not 
demonstrate distinct utilization patterns by Bacteroides 
species or host. The carbon utilization patterns had 
76.3% similarity across 190 carbon substrates, for all but 
six isolates (Fig. 4). At >90% similarity level, there were 
eight carbon utilization pattern groupings. Groups C2, 
C4, and C5 contained all three Bacteroides species from at 
least three hosts, groups C3, C6, and C7 contained B. 
thetaiotaomicron and B. xylanisolvens from at least two 
hosts, group C8 contained B. ovatus and B. xylanisolvens 
from all four hosts and finally grouping CI contained B 
thetaiotaomicron from human enrichments. Groups C4 
and C5 contained the most isolates with 31 in each 
(C4: 1, 16, and 14 and C5: 6, 1, and 24 numbers of 
B. ovatus, B. thetaiotaomicron, and B. xylanisolvens, 
respectively) The type strains B. thetaiotaomicron VPI 
5482 and B. xylanisolvens XBIA are found in the same 
group (C2) and had 93.7% similarity for carbon utiliza- 
tion patterns. The type strain B. ovatus ATCC 8483 did 
not group with any other isolates but was 82.2% similar 
to the other type strains. Two isolates (NLAE-zl-G481 
and NLAE-zl-P32) have low percent similarity to the 
other isolates, 41.9% and 30.3%, respectively, due to posi- 
tives across all 95 carbon substrate on a plate; NLAE-zl- 
G481 on plate PM2 and NLAE-zl-P32 on plate PMl. All 
of the Bacteroides isolates were able to utilize 11 carbon 
substrates: galactose, mannose, xylose, psicose, lyxose, 
allose, arabinose, 2-deoxy-D-ribose, 3-methylglucose, 
tagatose, and 5-keto-D-gluconic acid (Tables S2, S3). 
Bacteriodes ovatus isolates had the broadest range of car- 
bon substrate utilization with an additional 20 substrates 
utilized by all strains: fructose, a-D-glucose, maltose, 
melibiose, a-D-lactose, lactulose, sucrose, D-fructose-6- 
phosphate, maltotriose, glucuronamide, n-acetyl-D-galac- 
tosamine, 3-0-/J-D-galactopyranosyl-D-arabinose, and 
/J-methyl-D-galactoside. Bacteriodes thetaiotaomicron and 
B. xylanisolvens were all able to utilize three (rhamnose, 
maltotriose, and glucosamine) and four (rhamnose, fruc- 
tose, a-D-glucose, and glucuronamide) additional carbon 
substrates, respectively. 
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Figure 2. Cluster analysis of BOX primer rep-PCR. Similarity values (%). (A) Bacteroides ovatus; (B) S. thetaiotaomicron. Differentiation among 
host by isolate name: cow = NLAE-zl-Cx, goat = NI_AE-zl-Gx, Inuman = NLAE-zl-Hx, and pig = NLAE-zl-Px. 



230 



) 2014 The Authors. MicrobiologyOpen published by John Wiley & Sons Ltd. 



T. Atherly & C. J. Ziemer 



Substrate Utilization of Bacteroides Isolates 
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Figure 3. Cluster analysis of ERIC primers rep-PCR. Similarity values (%). (A) Bacteroides ovatus; (B) 6. thetaiotaomicron; (C) 6. xylanisolvens. Bar 
at top is similarity continuum. Differentiation among host by isolate name: cow = NLAE-zl-Cx, goat = NLAE-zl-Gx, human = NLAE-zl-Hx, and 
pig = NLAE-zl-Px. 
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Figure 3b. (Continued) 



232 



) 2014 The Authors. MicrobiologyOpen published by John Wiley & Sons Ltd. 



T. Atherly & C. J. Ziemer 



Substrate Utilization of Bacteroides Isolates 



NLAE-ZI-G37 
NLAE-ZI-H23 



xylanisolvens 
thetaiotaomicron 




^<C] B. thetaiotaomicron: NLAE-zl-H463, NLAE-zl-H492 
B. ovotus: NLAE-zl-HlSB, NLAE-zl-H251 
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B. thetaiotaomicron: NLAE-zl-C15, NLAE-zl-H207 
B. xylanisoivens: NLAE-zl-C233 

B. ovatus: NLAE-zl-H73 
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C504, NLAE-ZI-H409, NLAE-zl-P696, NLAE-zl-P737, NLAE-zl- 
C557, NLAE-ZI-H39, NLAE-zl-C484, NLAE-zl-C425, NLAE-zl- 
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H536, NLAE-ZI-H353 

B. xylanisoivens: NLAE-zl-C29, NLAE-zl-H194, NLAE-zl-P218, 
NLAE-ZI-P295, NLAE-zl-P402, NLAE-zl-P349, NLAE-zl-P887, 
NLAE-ZI-P352, NLAE-zl-P718, NLAE-zl-P371, NLAE-zl-P732, 
NLAE-ZI-P299, NLAE-zl-P393, NLAE-zl-P412 
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NLAE-zl-GlOg, NLAE-ZI-G408, NLAE-zl-G275, NLAE-zl-G353, 
NLAE-ZI-G421, NLAE-zl-G157, NLAE-zl-G44, NLAE-zl-G193, 
NLAE-ZI-G2, NLAE-ZI-G200, IMLAE-zl-GllO, NLAE-zl-G435, 
NLAE-ZI-G312 
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B. xylanisoivens: NLAE-zl-G39, NLAE-zl-P736 

B. ovatus: NLAE-zl-H356, NLAE-zl-H426, NLAE-zl-H313, NLAE- 
zl-Cll, NLAE-ZI-H255, NLAE-zl-H59 

B. xylanisoivens: NLAE-zl-H40, NLAE-zl-P225, NLAE-zl-G310, 
NLAE-ZI-H465 
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Figure 4. Cluster analysis of carbon substrate plates PM1 and PIVI2A utilization. Similarity values (%). Bar at top is similarity continuum. 
Differentiation among Inost by isolate name: cow = NLAE-zl-Cx, goat = NLAE-zl-Gx, human = NLAE-zl-Hx, and pig = NLAE-zl-Px. 
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S. ovatus: NLAE-zl-Cll, NLAE-zl-C34, NLAE-zl-C501, NLAE-zl-H153, NLAE-zl-H251, NLAE-zl- 
H304, NLAE-zl-HBSl, NLAE-zl-H541, NLAE-zl-H545, NLAE-zl-H59 

B. thetaiotaomicron. NLAE-zl-C13, NLAE-zl-C425, NLAE-zl-C484, NLAE-ZI-C504, NLAE-zl- 
C516, NLAE-ZI-C523, NLAE-zl-C557, NLAE-zl-G288, NLAE-zl-G295, NLAE-zl-G303, NLAE-zl- 
G493, NLAE-ZI-H207, NLAE-zl-H23, NLAE-zl-H334, NLAE-zl-H373, NLAE-zl-H39, NLAE-zl- 
H409, NLAE-ZI-H421, NLAE-zl-H486, NLAE-zl-H536, NLAE-zl-P696, NLAE-zl-P699, NLAE-zl- 
P737 

B. xylanisolvens: NLAE-zl-C139, NLAE-zl-C182, NLAE-zl-C195, NLAE-zl-C29, NLAE-zl-C315, 
NLAE-ZI-C339, NLAE-zl-G109, NLAE-zl-GllO, NLAE-zl-G157, NLAE-zl-G193, NLAE-zl-G194, 
NLAE-ZI-G2, NLAE-ZI-G200, NLAE-zl-G237, NLAE-zl-G275, NLAE-zl-G310, NLAE-zl-G312, 
NLAE-ZI-G346, NLAE-zl-G353, NLAE-zl-G37, NLAE-zl-G39, NLAE-zl-G397, NLAE-zl-G445, 



Nl 



NLAE-ZI-P352 xylanisolvens 

B. ovatus. NLAE-ZI-H255, NLAE-zl-H73 

NLAE-ZI-P732 xylanisolvens 

NLAE-ZI-P32 thetaiotaomicron 

^2:2-0 B. thetaiotaomicron. NLAE-zl-P750; B. xylanisolvens. NLAE-zl-P371 

B xylanisolvens XB1A xylanisolvens 

NLAE-ZI-P887 xylanisolvens 



N2 



N3 




B. ovatus. NLAE-ZI-H307, NLAE-zl-H313, NLAE-zl-H365, NLAE-zl-H426 

B. thetaiotaomicron. NLAE-zl-C15, NLAE-zl-G234, NLAE-zl-H353, NLAE-zl-H492, NLAE- 
ZI-P174, NLAE-ZI-P745 N4 

B. xylanisolvens. NLAE-zl-C178, NLAE-zl-C233, NLAE-zl-C367, NLAE-zl-G228, NLAE-zl- 
G289, NLAE-ZI-G405, NLAE-zl-G421, NLAE-zl-G44, NLAE-zl-P412 

B. ovatus. NLAE-ZI-H422; B. thetaiotaomicron: NLAE-zl-H463; B. xylanisolvens. NLAE-zl- |\|5 
C257 

NLAE-ZI-P299 xylanisolvens 
<] B. ovatus. B. ovatus ATCC 8483; 6. thetaiotaomicron: B. thetaiotaomicron VPI 5482 |\|6 
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NLAE-ZI-G435 xylanisolvens 



Figure 5. Cluster analysis of nitrogen substrate plate PM3B utilization. Similarity values (%). Bar at top is similarity continuum. Differentiation 
among host by isolate name: cow = NLAE-zl-Cx, goat = NLAE-zl-Gx, human = NLAE-zl-Hx, and pig = NLAE-zl-Px. 
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Using >90% similarity to assess nitrogen utilization 
patterns across 95 nitrogen substrates resulted in six 
groupings (Fig. 5). All three Bacteroides species from the 
four hosts were found in Nl, N4, and N5 groupings, 
group N3 contained B. thetaiotaomicron and B. xylanisol- 
vens from pigs, N2 contained B. ovatus from humans and 
N6 contained the B. ovatus and B. thetaiotaomicron type 
strains at 94.5% utilization similarity. Group Nl contained 
the largest number of isolates with 10 B. ovatus, 23 B. 
thetaiotaomicron, and 23 B. xylanisolvens. Isolate NLAE-zl- 
G435 had an extremely low similarity to the others, 4.5%, 
due to positives across all 95 nitrogen substrates. There 
were no nitrogen substrates that could be utilized by all 
of these Bacteroides isolates (Table S4). The nitrogen sub- 
strates most utilized were the purine derivatives xanthine, 
a nucleobase, and guanine. All B. ovatus and B. thetaiota- 
omicron isolates utilized xanthine, all B. ovatus isolates 
also utilized guanine; however, there were no nitrogen 
substrates that were utilized by all B. xylanisolvens isolates. 
The amino sugars D-galactosamine and D-mannosamine 
(B. ovatus 5 and 6, B. thetaiotaomicron 9 and 12, B. xylan- 
isolvens 15 and 14 isolates, respectively), N-acetyl-D-glu- 
cosamine and N-acetyl-D-galactosamine (B. ovatus 5 and 
5, B. thetaiotaomicron 14 and 10, B. xylanisolvens 19 and 
21 isolates, respectively), as well as the amino acids L-ho- 
moserine and L-ornithine (B. ovatus 4 and 4, B. thetaiota- 
omicron 9 and 7, B. xylanisolvens 12 and 10 isolates, 
respectively), were used by a moderate number of our 
isolates. 

Discussion 

Most of the research involving Bacteroides species has 
been done on isolates taken from humans (Salyers et al. 
1977; Salyers 1984; Xu et al. 2007; Zocco et al. 2007; 
Martens et al. 2009). In this study, we wanted to compare 
Bacteroides isolates within the B. thetaiotaomicron-ovatus- 
xylanisolvens clade obtained from different mammalian 
hosts. Our isolation of 103 bacteria within this clade 
demonstrated their presence in the GIT of mammals with 
varying diets and digestive physiologies. However, the 
Bacteroides species in this clade were not equally distrib- 
uted among host species. Bacteroides ovatus was recovered 
primarily in human and cow fecal enrichments. The 
B. xylanisolvens were isolated from all host fecal enrich- 
ments, but they were the most abundant from the goat 
fecal enrichments and least abundant from the human 
fecal enrichments. The B. thetaiotaomicron was also iso- 
lated from all host fecal enrichments with greatest abun- 
dance in isolates from pig fecal enrichments and relatively 
even distribution among other hosts. The primary host 
difference was more B. ovatus and B. thetaiotaomicron 
and less B. xylanisolvens isolates were from humans while 



in isolates from cow, goat, and pig B. xylanisolvens iso- 
lates predominated. 

Host-specific differences were expected for rep-PCR 
patterns as well as carbon and nitrogen substrate utiliza- 
tion. Based on analysis of fecal microbiota in 60 mam- 
mals using community 16S rRNA gene sequences 
indicating the influence of host diet type and digestive 
physiology (Ley et al. 2008a,b), there was potential to see 
host-specific or ruminant verses nonruminant patterns in 
the 16S rRNA gene sequences of our Bacteroides isolates. 
Only rep-PCR patterns supported genomic differences 
among B. ovatus, B. thetaiotaomicron, and B. xylanisolvens 
isolates due to host origin, indicative of coevolution in 
intestinal bacteria with their host. Neither the 16S rRNA 
gene sequences nor carbon and nitrogen substrate utiliza- 
tion were affected by the host they were isolated from. 
That 16S rRNA gene sequences do not reflect differences 
in host origin is not surprising due to its inherent proper- 
ties: universal distribution and function as well as high 
sequence conservation (Woese 2006). The specific envi- 
ronment and other members of the resident microbial 
community would influence mutations and horizontal 
gene transfer that could lead to genomic differences 
across strains of the same bacterial species. The rep-PCR 
data demonstrate host-specific patterns within the ge- 
nomes of the B. ovatus-thetaiotaomicron-xylanisolvens 
clade but the phenotypic array data suggest that the dif- 
ferences are not present in the carbohydrate utilization 
genes or if differences are there they maintain functional 
fidelity. 

While phylogentic grouping due to host origin was 
evident for B. ovatus isolates, there was no host origin 
distinction among B. thetaiotaomicron and B. xylanisol- 
vens 16S rRNA gene sequences. The B. thetaiotaomicron 
and B. xylanisolvens were obtained from all of the host 
fecal enrichments, while B. ovatus was obtained only 
from human and cow. Bacteroides ovatus isolates from 
cow and human had 2.0-2.8% difference in 16S rRNA 
gene sequence similarity, approaching the difference to 
separate bacterial species. The two human B. ovatus iso- 
lates, NLAE-Z1-H59 and NLAE-zl-H73 were only 96.5% 
similarity to the other B. ovatus and may represent a 
different Bacteroides species. In the context of this study, 
the 16S rRNA gene sequence was not able to determine 
the host origin of B. thetaiotaomicron and B. xylanisol- 
vans isolates. 

The rep-PCR method targets repetitive DNA sequences 
that are interspersed in the genome of the organism 
(Versalovic et al. 1994) and is used to visualize difference 
in the genomes. The amplification of these products can 
then be used to differentiate among strains of bacterial 
species. The Bacteroides are known to have mobile 
genetic elements such as conjugative transposons. 
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plasmids, and mobilizable transposons (Salyers 1984; 
Wexler 2007). These elements can allow for the uptake of 
foreign DNA from other microbes in the environment 
and possibly from the host (Comstock and Coyne 2003). 
On the basis of this information and reports of microbial 
communities across multiple host types (Dethlefsen et al. 
2007; Ley et al. 2008a), we anticipated differences in rep- 
PCR banding patterns due to mobilizable genetic ele- 
ments and host origin. Repetitive sequence-based PGR 
banding patterns grouped B. ovatus and B. xylanisolvens 
isolates by host origin. The B. thetaiotaomicron banding 
patterns resulted in groups of isolates from single and 
multiple host origin. The evidence of rep-PCR banding 
patterns grouping by host origin gives credence to a 
coevolution of the bacteria and host theory for B. ovatus 
and B. xylanisolvens, although less so for B. thetaiotaomi- 
cron. 

PhenotypicaUy, isolates from the three Bacteroides spe- 
cies shared the ability to utilize a number of the same 
carbon substrates, a reflection of their broad carbohydrate 
utilization. These results demonstrate functional redun- 
dancy; a predominant feature of the GIT microbiota (Ley 
et al. 2006), which would benefit the host by maintaining 
functionality of the intestinal microbiota even with the 
loss of a microbial lineage. This functional redundancy is 
considered to be host-driven selection (Ley et al. 2006), 
however, within the B. thetaiotaoniicron-ovatus— xylanisol- 
vens clade, host origin did not account for specific carbon 
utilization patterns. We speculate that the carbohydrate 
utilization functions of these Bacteroides are so important 
for utilization of plant cell wall carbohydrates that they 
are conserved across all mammalian intestinal microbiota. 
Taken together with the rep-PCR data, we also propose 
that the host-specific differences in the genomes of these 
Bacteroides species primarily come from genes not 
involved in carbohydrate utilization or any differences in 
carbohydrate utilization genes maintain precise function- 
ality. 

Our Bacteroides isolates showed an inability to utilize 
most nitrogen substrates on the nitrogen phenotypic 
array. Our results are supported by earlier findings that 
B. fragilis were unable to utilize free amino acids, 
peptides, nitrate, or urea as a nitrogen source (Varel and 
Bryant 1974). The lack of ammonia utilization by these 
isolates on the Biolog nitrogen phenotypic array was 
unexpected; thus, we used the media of Varel and Bryant 
(1974) to verify the results. Initial results demonstrated 
that the isolates could grow on the ammonia containing 
medium but also on the medium without ammonia; they 
were able to utilize the cysteine in the medium reducing 
agent (cysteine-HCL). This indicates that for results in 
the Varel and Bryant (1974) nitrogen source are con- 
founded with the cysteine in the reducing agent. We then 



repeated the test using sodium-thioglycolate (Sigma- 
Aldrich) as the reducing agent which confirmed that iso- 
lates could utilize ammonia and cysteine. This makes two 
nitrogen substrates that gave false negatives on the Biolog 
nitrogen phenotypic array. While we did not have con- 
cerns about false-positive reactions, there were no abiotic 
reactions for any nitrogen substrates, we utilized the 
sodium-thioglycolate-modified medium of Varel and 
Bryant (1974) to examine growth on animo sugars and 
purine derivatives using a subset of isolates and the type 
strains. The isolates were able to grow on D-mannos- 
amine hydrochloride, N-acetyl-D-glucosamine, N-acetyl- 
D-glactosamine, and xanthine, guanine but failed to grow 
on L-ornithine monohydrochloride (a potential false posi- 
tive on the Biolog plate). This supports the results of uti- 
lization of purine derivatives and amino sugars for some 
of the bacteria in the B. thetaiotaomicron-ovatus— xylan- 
isolvens clade. Macfarlane and Gibson (1991) reported 
that B. fragilis was able to grow on the glycoprotein 
mucin as well as mucin components such as N- acetylglu- 
cosamine and N-acetylgalactosamine. We have sequenced 
the genomes of 24 isolates and they all contained one 
copy of glutamine synthase and 2-3 copies of glutamate 
dehydrogenase (T. Atherly and C. J. Ziemer, unpubl. 
data), which appear to be functional based on our in 
vitro utilization assays with the amino sugars and purine 
derivatives. These results indicate that the Biolog PM3 
nitrogen phenotypic array is not a suitable method for 
determining nitrogen utilization of bacteria under anaero- 
bic conditions. Nitrogen utilization in the Bacteroides has 
not been well-defined, the Varel and Bryant (1974) paper 
is the primary source, and most research has focused on 
B. fragilis. A more detailed analysis of nitrogen substrate 
utilization by the Bacteroides, using an alternative reduc- 
ing agent to cysteine-HCL, is warranted. 

Isolates from the B. thetaiotaomicron-ovatus-xylanisol- 
vens clade were found in fecal enrichments fi-om mam- 
mals of differing digestive physiology, both ruminant and 
monogastric. Specificity of a particular Bacteroides species 
to one host, gut physiology or diet was not seen. The 
rep-PCR analysis indicated that the isolates' genomes had 
been altered by the host environment to some degree but 
genomic differences do not appear to be in 16S rRNA or 
carbohydrate utilization genes. These Bacteroides isolates 
were able to utilize a wide range of the carbon substrates 
but very few of the nitrogen substrates examined. Plant 
cell wall carbohydrates, recalcitrant to mammalian diges- 
tive enzymes, are the primary carbon substrates in the 
hind gut of both ruminants and monogastrics. The lack 
of host-specific carbon and nitrogen utilization in our 
Bacteroides isolates indicates conservation of functionality 
across hosts. This is the first report to indicate that host- 
associated genomic differences are outside of 16S rRNA 
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and carbohydrate utilization genes and to suggest conser- 
vation of specific bacterial species with the same fiinc- 
tionality across mammalian hosts for this Bacteroidetes 
clade. 

Acknowledgments 

This research was supported by a grant from Defense 
Advanced Research Projects Agency as part of it Intestinal 
Fortitude Program to C. J. Ziemer. The authors would 
like to thank Kerrie Franzen, John Tenhundfeld, and Ste- 
phen White for bacterial isolations. 

Conflict of Interest 

None declared. 
References 

Backhed, F., R. E. Ley, J. L. Sonnenburg, D. A. Peterson, and 
J. I. Gordon. 2005. Host-bacterial mutualism in the human 
intestine. Science (New York, N. Y.) 307:1915-1920. 

Bergman, E. N. 1990. Energy contributions of volatile fatty 
acids from the gastrointestinal tract in various species. 
Physiol. Rev. 70:567-590. 

Borglin, S., D. Joyner, J. Jacobsen, A. Mukhopadhyay, and T. 
C. Hazen. 2009. Overcoming the anaerobic hurdle in 
phenotypic microarrays: generation and visualization of 
growth curve data for Desulfovibrio vulgaris Hildenborough. 
J. Microbiol. Methods 76:159-168. 

Cole, J. R., Q. Wang, E. Cardenas, J. Fish, B. Chai, R. J. Farris, 
et al. 2009. The Ribosomal Database Project: improved 
alignments and new tools for rRNA analysis. Nucleic Acids 
Res. 37:141-145. 

Comstock, L. E., and M. J. Coyne. 2003. Bacteroides 
thetaiotaomicron: a dynamic, niche-adapted human 
symbiont. BioEssays 25:926-929. 

Dethlefsen, L., M. McFall-Ngai, and D. A. Relman. 2007. An 
ecological and evolutionary perspective on human-microbe 
mutualism and disease. Nature 449:811-818. 

Eckburg, P. B., E. M. Bik, C. N. Bernstein, E. Purdom, M. 
Sargent, S. R. Gill, et al. 2006. Diversity of the human 
intestinal microbial flora. Science 308:1635-1638. 

Flint, H. J., E. A. Bayer, M. T. Rincon, R. Lamed, and B. A. 
White. 2008. Polysaccharide utilization by gut bacteria: 
potential for new insights from genomic analysis. Nat. Rev. 
Microbiol. 6:121-131. 

Lane, D. J. 1991. 16S/23S rRNA sequencing. Pp. 115-163 in 
E. Stackebrandt and M. GoodfeUow, eds. Nucleic acid 
techniques in bacterial systematics. Willey, New York, NY. 

Ley, R. E., D. A. Peterson, and J. I. Gordon. 2006. Ecological 
and evolutionary forces shaping microbial diversity in the 
human intestine. Cell 124:837-848. 



Ley, R. E., M. Hamady, C. Lozupone, P. J. Turnbaugh, R. R. 

Ramey, J. S. Bircher, et al. 2008a. Evolution of mammals 

and their gut microbes. Science 320:1647-1652. 
Ley, R. E., C. A. Lozupone, M. Hamady, R. Knight, and J. I. 

Gordon. 2008b. Worlds within worlds: evolution of the 

vertebrate gut microbiota. Nat. Rev. Microbiol. 6:776-788. 
Macfarlane, G. T., and G. R. Gibson. 1991. Formation of 

glycoprotein degrading enzymes by Bacteroides fragilis. 

FEMS Microbiol. Lett. 61:289-293. 
Mappley, L. J., M. L. Black, M. AbuOun, A. C. Darby, M. J. 

Woodward, J. Parkhill, et al. 2012. Comparative genomics 

of Brachyspira pilosicoli strains: genome rearrangements, 

reductions and correlation of genetic compliment with 

phenotypic diversity. BMC Genomics 13:454, http://www. 

biomedcentral.com/ 1 47 1 -2 1 64/ 1 3/454. 
Martens, E. C, N. M. Koropatkin, T. J. Smith, and J. I. 

Gordon. 2009. Complex glycan catabolism by the human 

gut microbiota: the Bacteroidetes Sus-like paradigm. J. Biol. 

Chem. 284:24673-24677. 
Martens, E. C, E. C. Lowe, H. Chiang, N. A. Pudlo, M. Wu, 

N. P. McNulty, et al. 2011. Recognition and degradation of 

plant cell wall polysaccharides by two human gut symbionts. 

PLoS Biol 9:el001221. 
Oh, P. L., A. K. Benson, D. A. Peterson, P. B. Patil, E. N. 

Moriyama, S. Roos, et al. 2010. Diversification of the gut 

symbiont Lactobacillus reuteri as a result of host-driven 

evolution. ISME J. 4:377-387. 
Salyers, A. A. 1984. Bacteroides of the human lower intestinal 

tract. Annu. Rev. Microbiol. 38:293-313. 
Salyers, A. A., J. R. Vercellotti, S. E. West, and T. D. Wilkins. 

1977. Fermentation of mucin and plant polysaccharides by 

strains of Bacteroides from the human colon. Appl. Environ. 

Microbiol. 33:319-322. 
Thomas, F., J. H. Hehemann, E. Rebuffet, M. Czjzek, and 

G. Michel. 2011. Environmental and gut bacteroidetes: the 

food connection. Front. Microbiol. 2:1-10. 
Varel, V. H., and M. P. Bryant. 1974. Nutritional features of 

Bacteroides fragilis suhsp. fragilis. Appl. Microbiol. 18:251- 

257. 

Versalovic, J., T. Koeuth, and J. R. Lupski. 1991. Distribution 
of repetitive DNA sequences in eubacteria and application 
to fingerprinting of bacterial genomes. Nucleic Acids Res. 
19:6823-6831. 

Versalovic, J., M. Schneider, F. J. de Bruijn, and J. R. Lupski. 

1994. Genomic fingerprinting of bacteria using repetitive 

sequence-based polymerase chain reaction. Methods Mol. 

Cell. Biol. 5:25-40. 
Wexler, H. M. 2007. Bacteroides: the good, the bad, and the 

nitty-gritty. Clin. Microbiol. Rev. 20:593-621. 
Woese, C. R. 2006. A new biology for a new century. 

Microbiol Mol. Biol. Rev. 68:173-186. 
Xu, J., M. A. Mahowald, R. E. Ley, C. A. Lozupone, M. 

Hamady, E. C. Martens, et al. 2007. Evolution of symbiotic 

bacteria in the distal human intestine. PLoS Biol. 5:el56. 



2014 The Authors. MicrobiologyOpen published by John Wiley & Sons Ltd 



237 



Substrate Utilization of Bacteroides Isolates 

Zaneveld, J., P. J. Turnbaugh, C. Lozupone, R. E. Ley, J. I. 

Gordon, and R. Knight. 2009. Host-bacterial coevolution 

and the search for new drug targets. Curr. Opin. Chem. 

Biol. 12:109-114. 
Ziemer, C. J. 2013. Broad diversity and newly cultured 

bacterial isolates from enrichment of pig feces on complex 

polysaccharides. Microb. Ecol. 66:448-461. 
Zocco, M. A., M. E. Ainora, G. Gasbarrini, and A. Gasbarrini. 

2007. Bacteroides thetaiotaomicron in the gut: molecular 

aspects of their interaction. Dig. Liver Dis. 39:707-712. 

Supporting Information 

Additional Supporting Information may be found in the 
online version of this article: 



T. Atherly & C. J. Ziemer 

Table SI. Bacteroides isolates with NCBI Accession Num- 
bers by host origin and Bacteroides species. 
Table S2. List of phenotypic microarray substrates and 
the number of isolates by Bacteroides species that utilized 
the substrates of Biolog PMl. 

Table S3. List of phenotypic microarray substrates and 
the number of isolates by Bacteroides species that utilized 
the substrates of Biolog PM2A. 

Table S4. List of phenotypic microarray substrates and 
the number of isolates by Bacteroides species that utilized 
the substrates of Biolog PM3B. 



238 



2014 The Authors. MicrobiologyOpen published by John Wiley & Sons Ltd. 



